PITHA 07/03 
TTP/07-11 
SFB/CPP-07-22 
18 June 2007 



Ultrasoft contribution to quarkonium production 

and annihilation 



M. BENEKE a , Y. KlYO 6 , A. A. PENIN 6 ' c ' d 

a Institut fur Theoretische Physik E, RWTH Aachen, 
D-52056 Aachen, Germany 

b Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, 
D-76128 Karlsruhe, Germany 

c Department of Physics, University Of Alberta 
Edmonton, AB T6G 2J1, Canada 
d Institute for Nuclear Research, Russian Academy of Sciences, 
119899 Moscow, Russia 



Abstract 

We compute the third-order correction to electromagnetic S-wave quarkonium 
production and annihilation rates due to the emission and absorption of an ul- 
trasoft gluon. Our result completes the analysis of the non-relativistic quarkonium 
bound-state dynamics in the next-to-next-to-next-to-leading order. The impact of 
the ultrasoft correction on the T(1S) leptonic width and the top quark-antiquark 
threshold production cross section is estimated. 



PACS numbers: 12.38.Bx, 14.40. Gx, 14.65.Fy, 14.65.Ha 



1 Introduction 



The theoretical study of non-relativistic heavy quark-antiquark systems is among the 
earliest applications of perturbative quantum chromodynamics (QCD) pQ. Perturbation 
theory applies to the bound-state dynamics of bottomonium, at least within the sum rule 
approach [2] , and top-antitop systems [3] , since non-perturbative effects are under control 
[Hl5] . This makes heavy quark-antiquark systems well suited to determine fundamental 
parameters of QCD, the strong coupling constant a s and the heavy-quark masses m. 

The binding energy of a quarkonium state and the value of its wave function at the 
origin - field-theoretically, the residues of two-point functions of local currents - are of 
primary phenomenological interest. The former determines the mass of the bound state, 
while the latter controls its production and annihilation rates. The quarkonium ground- 
state energy has been computed through 0{ma 5 s ) including the third-order correction 
to the Coulomb approximation [BITT] . This result has been extended to the excited S- 
wave states jSE]- For the wave function at the origin a complete result is only available 
including the second-order correction [T0~Jll||12| . This correction is large even for top 
quarks, and for a reliable perturbative prediction the third-order approximation seems to 
be needed. This amounts to a difficult calculation, which can be broken into several well- 
defined pieces, some of which are already available, such as the double-logarithmically 
enhanced C(a^ln 2 a s ) terms [T31[T3] and the single-logarithmic 0(a^lna s ) terms [TofTo] . 
The calculation of the most difficult non-logarithmic term has been started in [819"] . 
where the contribution to the wave function at the origin from the loop corrections to 
the colour-Coulomb potential have been evaluated. In this and the companion paper [17] 
the remaining contributions from the non-Coulomb potentials and due to the emission 
and absorption of an ultrasoft gluon by the quarkonium bound state are presented. 
The ultrasoft correction discussed below is of special interest, because it constitutes a 
qualitatively new effect, which shows up for the first time in the third order. No other 
such effects are expected in higher orders of the perturbative expansion. The complete 
third-order correction to the wave function at the origin can now be expressed in terms of 
a few yet unknown matching coefficients, which can be obtained by standard fixed-order 
loop calculations. 

2 Ultrasoft correction to the wave function 
2.1 Definitions 

In non-relativistic bound states the quark velocity v is a small parameter. An expansion 
in v may be performed directly in the QCD Lagrangian by using the framework of 
effective field theory [18.19.20j, or diagrammatically with the threshold expansion [21 J . 
The relevant momentum regions are the hard region (energy k° and momentum k of 
order m), the soft region (k°,k ~ mv), the potential region (k° ~ mv 2 , k ~ mv), and 
the ultrasoft region (k°, k ~ mv 2 ). Integrating out the hard modes amounts to matching 
onto non-relativistic QCD (NRQCD) [19]. If one also integrates out the soft modes and 
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potential gluons, one obtains the effective theory called potential NRQCD (PNRQCD), 
which contains potential heavy quarks and ultrasoft gluons as dynamical fields [20J (see 
also [22J ) . In this theory the leading colour- Coulomb potential is part of the unperturbed 
Lagrangian, so that the propagation of a colour-singlet quark- ant iquark pair is described 
by the Green function of the Schrodinger equation 

(H -E)G<g\r,r';E) = SW(r-r'), (1) 

with 
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m 

r = \r\, m the heavy-quark pole mass, and C F = (iV c 2 -l)/(2iV c ), N c = 3. The PNRQCD 
Lagrangian further contains interactions of quarks with the multipole-expanded ultra- 
soft gluon field and instantaneous, spatially non-local interactions ("potentials"), which 
can be treated as perturbations. This constitutes the basic framework for the per- 
turbative analysis of quarkonium bound-state properties. The colour- singlet Coulomb 
Green function Gq (r, r'; E) has an infinite number of bound-state poles with energies 
E^ = —m(a s CF) 2 /(2n) 2 and contains the information about the corresponding wave 
functions. In the quark-antiquark Fock state sector the perturbations to the energy 
levels and wave functions can be taken into account by replacing Hq by the PNRQCD 
Hamiltonian H with the ultrasoft modes are excluded. 

In this paper, however, we are interested in the leading ultrasoft effect. To connect 
the concept of a non-relativistic wave function at the origin to a physical quantity, we 
consider the two-point function of the electromagnetic heavy-quark current = Qj^Q 
in full QCD, 

(qtf* ~ 9,uq 2 ) I% 2 ) = i J Ae v (0\T 3fl (x) Ju (0)\0) , (3) 

whose poles are related to electromagnetic production and annihilation rates of the 
corresponding quarkonium states. In PNRQCD j M is represented in terms of operators 
constructed from the non-relativistic quark and antiquark two-component Pauli spinor 
fields ip an d X? 

j = c^ax + -^^(tD 2 X + ■■■■ (4) 

The matching coefficients c„(/i) = 1 + c^(ju) a s (fi) / (An) + . . . and d v (/j) = 1 + 0(a s ) 
represent the contributions from the hard modes with /i a factorization scale that is also 
implicit in the renormalization convention for the operators on the right-hand side. (We 
use dimensional regularization with d = 4 — 2e and the MS scheme.) We also introduce 
the PNRQCD two-point function 

2(d - l)N c G(E) = i [ d d xe iE *° (0|T[^V X ](x) frVV](0)|0) , (5) 
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Figure 1: The ultrasoft correction as a PNRQCD Feynman diagram. The bold and the double 
lines stand for the singlet and octet Coulomb Green functions, respectively, the curly line 
represents the ultrasoft-gluon propagator, the black circles represent the chromoelectric dipole 
interaction g s rE, and the squares correspond to the non-relativistic currents. 



where E = ^fq 1 — 2m. In leading order G(E) coincides with (the correspondingly 
regularized) Gq\o, 0; E). Substituting the expansion (@J into ([3]) and using an equation- 
of-motion relation for the insertion of the derivative current in (HI), we obtain 



N r 
2m 2 




+ ... 



G{E) , (6) 



which is valid up to the third order. G(E) has Coulomb bound-state poles at energies 
E n sa with spin and orbital angular momentum 5 = 1 and I = 0, respectively, 
following from the form of the current, ip^crx- Near the pole 

G(0 A E)^ J M °f , (7) 

v ' ' ; E^E n E n - E-ie w 

which defines the "wave function at the origin" . Beginning from the second-order correc- 
tion, these wave-functions are factorization-scale dependent, but this dependence cancels 
in the residues of the poles of n(g 2 ), which determine the observable electromagnetic pro- 
duction and annihilation rates. 

The third-order (NNNLO) corrections to |^ n (0)| 2 originate (i) from single insertions 
of the third-order potentials in the PNRQCD Lagrangian and multiple insertions of flrst- 
and second-order potentials, where the order is determined by the combined suppression 
in a s and v relative to the leading-order Coulomb potential a s /r ~ a s v . This contribu- 
tion has been computed in [8,9,17] and expressed in terms of the yet unknown three- loop 
correction to the Coulomb potential and four constants related to (9(e) parts of loop cor- 
rections to other potentials; (ii) from the emission and absorption of an ultrasoft gluon 
(similar to the Lamb shift for energy levels), calculated below. 



2.2 Calculation of the ultrasoft correction 

The leading ultrasoft interactions in the PNRQCD Lagrangian are g s ip lf (x)(A (t, 0) — 
xE(t,0))tp(x) together with a similar term for the antiquark field. The contribution 
from the A coupling cancels (or can be gauged away), leaving the chromoelectric dipole 
interaction, which results in a NNNLO correction [23,24,25] . The PNRQCD diagram 
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Figure 2: The NRQCD decomposition of the PNRQCD vertex. The dashed line corresponds 
to a potential Coulomb gluon. The black circles on the right-hand side correspond to the 
NRQCD vertices g s pA/m. The bold arrows correspond to the potential quark and antiquark 
propagators with any number of the potential gluon exchanges. 

representing this correction is shown in Fig. [TJ The corresponding correction to G(E) 
denned in ([5]) reads 



with the understanding that one picks up only the pole at k° = \k\ — ie in the gluon 
propagator. Here G^\r,r';E) denotes the colour-octet Coulomb Green function corre- 
sponding to the Hamiltonian ([2]) with the octet potential Vq (r) = — Cp)a s /r 
(Ca = N c ). Only the / = 1 partial wave of the octet Green function contributes to (|SJ), 
and since the octet potential is repulsive, G c (r, r ; E) does not have bound-state poles. 

Expression (JSJ) cannot be used in practice, because the ultrasoft correction is di- 
vergent. Its definition requires specifying a regulator and subtractions, which must be 
chosen to be consistent with the calculation of the potential insertions; the potentials 
themselves; and the hard matching coefficients. In the following we provide the necessary 
definitions, deferring a detailed discussion of many interesting technical aspects of the 
calculation to a later publication that will provide results not only for |-^ n (0)| 2 but also 
for the full correlation function G(E). To apply dimensional regularization we transform 
(jHJ) to momentum space. We also find it convenient to re-express the PNRQCD vertex in 
terms of the NRQCD vertices from which it is derived by equation-of-motion relations, 
see [25] and Fig. [2] for a graphical representation. The formulation of non-relativistic 
effective theory in dimensional regularization PI12II26II27] is very convenient, because 
it allows one to combine bound-state calculations with loop calculations of matching 
coefficients, which are technically feasible only in dimensional regularization. Further- 
more, when loop integrals are expanded in the sense of the threshold expansion [21], the 
matching of contributions from different regions is automatic. 

The integral over the three-momentum k of the ultrasoft gluon is ultraviolet (UV) 
divergent. The divergence is related to the factorization of the ultrasoft scale from the 
other scales, and cancels when all pieces of the calculation are added. The UV-divergent 
part of the ultrasoft integral has the form of a single insertion of a third-order potential 
and of a one-loop correction to the coefficient d v of the derivative current in (jlj). In 
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(8) 
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Figure 3: Example of a three-loop vertex correction with an ultrasoft exchange. The right- 
hand side illustrates the subtraction of the UV divergent part. The thin arrows correspond to 
free potential quark and antiquark propagators. 



fact, they cancel precisely infrared (IR) divergences in the calculation of these quantities 
[6||28] . We therefore define the ultrasoft correction by adding counterterms that cancel 
these ultrasoft subdivergences, 



S US G{E) 



d d-i £ d d-i £ > 
(2?r) d - 1 (2tt) 



+ 



d d-i p d d-i p 



-Sd 



jG^^p-^^U + dV,,] G ( £>(p',£';E) \, (9) 



'div 



e + e 2 ~ 



6m 2 



(2ir) d - 1 (2Tr) d 

where 5U follows from (jSj), and the counterterms read (q = p — p' 



6V, 



c.t. 



a s Ci 
6e 



rt a s , A(r a ^ na > 



C^-f + 4(Ci + 2C A C F ) 



C A \ ol 



+16 C F --f K + 



a s p 2 + p' 2 



2g 2 



div 



a s 16Cf 
Air e 



(10) 



(11) 



The counterterms added here are subtracted from the other parts of the calculation |17j . 
With ft 2 = e lE fi 2 / (4tt) subtracting poles in e corresponds to the MS subtraction. 

After adding this counterterm the potential loops and Coulomb Green functions can 
be evaluated in three dimensions unless the potential loop integrations are divergent. 
Such divergences occur in potential vertex subgraphs up to three loops, and they cor- 
respond to IR divergences in the three- loop correction to c v . (There are also over-all 
divergences in G(E), but they are irrelevant to the calculation of the bound-state pole 
and residue.) We separate this vertex sub divergence by adding and subtracting the 
three-loop vertex subdiagram at zero external momentum p, as shown in Fig. [3j The 
vertex UV divergence is logarithmic, and does not depend on the external momentum. 
Hence, it is isolated in first term on the right-hand side of the equation of Fig. [3j while 
the difference in the brackets is finite and can be computed in three dimensions. The 
first term factorizes into a three- loop diagram at p — 0, which has to be computed in 
d dimensions, and the leading-order expression for G(E). The divergent part cancels 
the contribution from c v , and the remainder can again be evaluated in three dimensions. 
Thus, we do not need the d — 1 dimensional Coulomb Green function, which is unknown. 
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n 


1 


2 


3 


4 


5 


6 




353.06 


256.62 


224.26 


206.88 


195.48 


187.16 



Table 1: Numerical result for the non- logarithmic part of the ultrasoft contribution, as defined 
in (|13p ; n denotes the principal quantum number. 



2.3 Result 

We write the perturbative expansion for the wave function at the origin as 

l^n(0)| 2 = |^(0) | 2 (l + S^ n + 5^ n + . 



(12) 



where the leading-order Coulomb wave function at the origin in three dimensions is given 
by ^^(O)! 2 = (ma s CV) 3 /(87rn 3 ), and 5^ip n stands for the mth order correction. For 
the ultrasoft part of the third-order correction we obtain 



a'' 



71 



I 1 

-C\C F + — C A C 2 F + -C% 



1 

-( 

6 



1 



-C 2 A Cp + -^CaCp + -C| 



e m 



+ 



+ 



+ 



+ 



ln2 N 



18 



6 



25 



— — I CiC F + — - - In 2 C A Cl + — - In 2 C 



16 



12 6 

8 



19 



18 



ZL* a Uf —LyA^F ~ 



In a s + 



6 6 3 



ln 2 ^ 



-C a Cf + —Ca^f + 77^F 



In a s In — 

m 



52 8 



c a + - o ln2 - + 6 - — - - ln2 - — F n C A C 



+ 



9 



10 



32 



3n 2 



ln a. 



+ 



C\ + 



11 5 

y + 3 



+ ^ln2 + ^ff n )ClC F + 



3 
2 



3rr 



20 

y 



|5 + — -61n2 + -H n K' 



ln^ + 5: 



??? 



+ ^ln2 + =f#„)C A C 2 



(13) 



where H n = ln — - + X)^=i i- The mos t difficult part of the calculation is the non- 
logarithmic term which we could compute only numerically. For the six lowest states 
its value is given in Tabled) Note that \ip^(0) | 2 in (fT2l is formally defined in d dimensions, 
but as explained above, we do not need the explicit <i-dimensional expression, because 
the 1/e pole terms in (fT3"j) cancel with pole terms from c 2 \ip^ (0)| 2 contained in (jBJ). 
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We verified that the logarithmic terms in f|T3|) when combined with those from the 
third-order potentials insertions [T7] agree with [15^1 . The sum of the 1/e poles in the 
complete third-order correction to |-?/> n (0)| 2 from ([TBI and [T7] combined contains a term 
proportional to the single insertion of the first-order Coulomb potential and a term that 
must cancel against the infrared pole in twice the three-loop correction to c v . This second 
term reads 



—CaC f + — Cp ) ^ + 



1 

24 



1 

36 




CaC f 



«-ln2 

72 



^C A C F T F n f + ^C F T F -^C F T F n f 



(14) 



where (3q = 11CU/3 — 4T F nj/3 is the one- loop QCD beta-function, T F = 1/2 and rif is 
the number of light-quark flavors. The n y-part of the three-loop coefficient cffl is known 
[29] . and we checked that the n/-dependent pole parts cancel as required. Eq. (fi~4l is 
consistent with the scale dependence of the hard matching coefficient given by Eqs. (8), 
(9) of [15], except for the rational part of the CaC f term in 7^ 3 \ which should be 
increased by 7/24 to agree with our result D The difference is due to the fact that by 
ijk a k spin-algebra in [15] is not completely (^-dimensional, hence the 
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using [a , a J 

result for the scale-dependence of given there does not correspond to a calculation 
in conventional dimensional regularization. 



3 Quarkonium phenomenology 

We briefly discuss the size of the ultrasoft correction for the two most relevant cases, the 
leptonic decay of the T(IS') and the threshold production of top quark- antiquark pairs 
in e + e _ annihilation. This discussion must necessarily be preliminary, since the ultrasoft 
correction alone is factorization scheme and scale dependent. For the quarkonium spin- 
triplet ground state, n = 1, we obtain from ( [TBI , omitting the 1/e poles, 

5 US ^ = a 3 s | - 18.71 In 2 a, + 52.03 In a, + 112.38 

In ^-6.55 In 2 H. (15) 
m mj 



+ 



23.52 In a, - 30.98 



The scale of the coupling a s is most naturally of order of the inverse Bohr radius ma s C F 
in two of the three powers of the overall factor a 3 , and of order of the ultrasoft scale ma 2 s 



1 In Eq. (7) of p~5] for the excited states the term 4Cp(l — l/n 2 )/3 is missed and the Polnn term 
should be multiplied by three. 

2 In addition, as noted in [25], the term — |/3o7ii 2 ^ in Eq. (8) of [IS] should read +/3o7« 2 ^ 
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in the third. However, any other scale choice is formally equivalent at this order. In the 
following we evaluate a s at fiB = mCj?a s (/i£), wherever it appears. The scale \i in the 
ln(/x/m) terms is related to scale-dependent potentials and hard matching coefficients. 
We vary /i/m between a s Cp (corresponding to the scale \ib) and 1 (hard scale). 



T(IS') leptonic width. Up to a normalization factor the T(15) leptonic decay width is 
given by the residue of n(g 2 ) at the ground state pole. The leading-order expression for 
the decay width follows from resulting in 1^° = 47riV c e 2 a 2 |'?/>f' (0)| 2 / (3m 2 ), where 
e q is the electric charge of the quark flavour q and a is the fine-structure constant. 
The non-perturbative contribution to the width is quite sizeable and out of control for 
higher resonances [30 31J, hence we consider only the case n = 1. Adopting a s = 0.30, 
which corresponds to hb ~ 2 GeV for the T(IS'), we obtain 8T\/T\ ~ 3.0 from the 
non-logarithmic correction 5™ s alone. Proceeding as described above for the logarithmic 
terms results in the estimate 

5T 1 w [0.61 - 1.93] T\°. (16) 

It therefore appears that the large non-logarithmic term leads to a large enhancement of 
the width. Whether or not perturbation theory is out of control (as may be suggested by 
the upper limit of the given range) can be decided only after combining all third-order 
terms. 



Top-quark production near threshold. For top quarks non-perturbative effects are neg- 
ligible, but its decay width T t smears out the Coulomb resonances below the thresh- 
old. The NNLO analysis of the cross section [52] shows that only the ground-state 
pole gives rise to a prominent resonance. Although the calculation of the normal- 
ized cross section R = a(e + e~ — > tiX)/a(e + e~ — > /U + /i~) requires the full Green 
function G(E), the height of the resonance can be estimated from the wave function 
at the origin of the would-be toponium ground state. In the leading-order approxi- 
mation R\° « 67rN c ef\ipi (0)\ 2 / (m 2 T t ). Adopting a s = 0.14, which corresponds to 
fiB ~ 32.5 GeV, we obtain 5R\/ R\ ~ 0.31 from the non-logarithmic correction 8™ s alone. 
Including an estimate of the logarithmic terms we find 

8Rt w [(-0.17) - (+0.13)] R\°. (17) 

Hence, despite the large quark mass, we may anticipate a sizeable third-order correction, 
unless there are cancellations. 



4 Summary 

We evaluated the third-order correction to electromagnetic quarkonium production and 
annihilation due to the emission and absorption of an ultrasoft gluon. Together with 
other contributions [8"f9"|13lll4|ll5|ll6|ll7|l29j already completed the problem of evaluating 
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the total 0{a Z g) corrections is now reduced to the calculation of four 0(e) terms in the 
NNNLO heavy quark- ant iquark potential [17], the three-loop colour-singlet Coulomb 
potential, and the three-loop vector current matching coefficient in the MS scheme in 
dimensional regularization. The previously unknown non- logarithmic ultrasoft contribu- 
tion is large and significantly increases the production and annihilation rates. It might 
limit the accuracy of the perturbative analysis of the quarkonium even for top quarks. 
We should however emphasize that a definite conclusion can only be drawn once the full 
NNNLO result is available. In this respect the sizable negative third-order correction 
from the perturbation potentials [5119717] should be mentioned. 
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